Vibrational properties of graphene fluoride and graphane 
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The vibrational properties of graphene fluoride and graphane are studied using ab initio calculations. We 
find that both sp 3 bonded derivatives of graphene have different phonon dispersion relations and phonon 
density of states as expected from the different masses associated with the attached atoms fluorine and 
hydrogen, respectively. These differences manifest themselves in the predicted temperature behavior of the 
constant- volume specific heat of both compounds. 



PACS numbers: 73.22.-f, 73.22.Pr, 71.15.Mb 

The chemical modification of graphene 1 - is being the 
focus of increasing interest^— Radicals such as oxygen, 
hydrogen, or fluorine atoms can be adsorbed on the sur- 
face of graphene. The adsorbed radicals can attach to 
the graphene layer in a random way, as is the case in 
graphene oxide (GO)^ or they can form ordered patterns. 
The later is expected to be the case for hydrogen^J^ and 
fluorine adsorbates . 11 ' 12 The two-dimensional materials 
that form in those cases^ have been named graphan o 3 i 10 
and graphene fluoride 4 (or fluorographene 6 ), respectively. 
Raman spectroscopy provides useful information, and it 
is being actively used, during the hydrogenation 3 - and 
fluorination 4 ^ process of graphene. Both graphane and 
graphene fluoride are expected to be wide band gap 
materials J£ As the energies of the lasers used were lower 
than the band gap, there was no signature of Raman ac- 
tivity from fully hydrogenated or fully fluorinated regions 
in those experiments and the Raman peaks were as- 
sociated with graphene. 

The main features in the Raman spectra of pris- 
tine graphene is the appearance of the G and 2D 
bands, which appear around 1580 cm -1 and 2680 cm -1 , 
respectivelyi^ 1 ^ The G-band is associated with the dou- 
bly degenerate E2 9 phonon mode of graphene, while the 
2D mode (also called G') originates from a second-order 
process, involving two phonons near the K point without 
the presence of any kind of disorder or defect^ On the 
other hand, the presence of defects in the sample (such 
as adatoms in partially hydrogenated or fluorinated sam- 
ples) activates additional peaks in the Raman spectra 
of graphene. This is the case of the D, D', and D+D' 
peaks that appear at 1350 cm -1 , 1620 cm -1 , and 2950 
cm -1 , respectively These Raman peaks originate from 
double-resonance processes at the K point in the presence 
of defects. As the D, D and D+D' peaks involve phonon 
modes from graphene, they appear at the same frequen- 
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cies in partially hydrogenated 3 or partially fluorinated 
samples £ 

The D and G peaks of graphene provide valuable infor- 
mation on the density of defects. These peaks tend to dis- 
appear in the limit of almost fully covered samples as was 
recently found in graphene fluoride&i but not in previous 
experiments on graphane^ or graphite fluorides i 14 ' 15 The 
first report on Raman signatures of multi-layer graphene 
regions that are fully covered with fluorine atoms was 
presented in Ref. By using a UV laser with an en- 
ergy of 5.08 eV, two Raman active modes were detected 
at 1270 cm -1 and 1343 cm -1 , which are absent for lower 
laser energies. An infrared active mode was also reported 
at 1204 cm -1 . These modes were correlated with DFT 
calculations of the phonon frequencies of graphene fluo- 
ride at the T point 4 

In this letter, by using first-principles calculations, we 
investigate the phonon dispersion and the phonon den- 
sity of states of graphene fluoride and graphane. This 
information can be useful for the interpretation of future 
experiments on infrared, Raman, and neutron-diffraction 
spectra of these novel two-dimensional compounds, as 
well as in the study of a wide variety of other physical 
properties such as specific heat, thermal expansion, heat 
conduction, and electron-phonon interaction. 
All reported results are obtained with the ABINIT 
code^ using a plane wave basis and Troullier-Martins 
pseudopotentialSi— An energy cutoff of 40 Ha was used. 
A 30 Bohr layer of vacuum was used to separate the 
graphane/fluorographene sheet from its periodical im- 
ages in order to avoid unphysical interactions. The Bril- 
louin Zone was sampled using a 30x30x1 Monkhorst-Pack 
gricUi for the electronic structure. Both considered sys- 
tems are fully relaxed such that all forces are smaller 
than 5 x 10 -5 Hartree/Bohr and all stresses smaller 
than 5 x 10~ 7 Hartree/Bohr 3 . The dynamical proper- 
ties are calculated within the density-functional pertur- 
bation theory 1 ^ (DFPT), as this allows us to calculate 
arbitrarily phonon frequencies without requiring the use 
of a supercell. 

Figure[T]shows the calculated phonon dispersion relation 
for the chair conformation of graphane, as this is the en- 
ergetically most favourable crystal configuration^ The 
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FIG. 1. (Color online) Phonon dispersion of graphane in 
the chair conformation. The dots are the directly calculated 
frequencies and the lines are interpolated values. The inset 
shows the first Brillouin zone and the wavevector path used. 
The right-hand side shows the phonon DOS, where the black 
lines indicate the total DOS, the blue dotted lines the C pro- 
jected DOS, and the red dashed line the H projected DOS. 



FIG. 2. (Color online) Phonon dispersion of fluorographene 
in the chair conformation. The dots indicate the directly cal- 
culated frequencies and the lines are the interpolated values. 
The right-hand side shows the phonon DOS, where the black 
lines indicate the total DOS, the blue dotted lines the C pro- 
jected DOS, and the green dashed line the F projected DOS. 



dots indicate the directly calculated frequencies, while 
the solid line is obtained by using a Fourier based inter- 
polation of the interatomic force constants obtained on a 
10x10x1 grid of ^-points. As can be seen, the agreement 
between both calculations is very good. The absence of 
any imaginary frequency indicates that the chair con- 
formation is stable. It is apparent that the phonons of 
graphane can be divided into low-, intermediate-, and 
high-frequency groups of phonons. The right-hand side 
of the figure shows the corresponding density of states 
(DOS), obtained with the interpolated frequency values. 
The blue dotted line is the DOS projected on the C 
atoms, and the red dashed line is the projection on the H 
atoms. From the projected phonon DOS, one can imme- 
diately identify the high-frequency modes as dominantly 
H modes, as can be expected from the C-H stretching 
modes. The acoustic modes are dominantly C modes 
while the intermediate-frequency group of phonons have 
significant contributions from both types of atoms. 

A similar approach can be applied to fluorographene, 
resulting in the phonon dispersion relation (left) and pro- 
jected phonon DOS (right) shown in Fig. [5J In the pro- 
jected DOS, the blue dotted line is the projected DOS on 
the C atoms, and the green dashed line is the projection 
on the F atoms. From the dispersion we found that also 
the chair conformation of fluorographene is stable. 
There arc substantial differences in the phonon spectra 
and DOS of graphane and fluorographene. The latter 
does not show clearly separated groups of phonons. The 
high-frequency modes are also lower in frequency in com- 



parison with similar modes in graphane because F is 
heavier than H. In this case, the dominant contribution 
to the acoustic modes comes from the F atoms while the 
high-frequency modes around 1200 cm" 1 have a clear C 
character. This behavior is the opposite from the one in 
graphane because the mass of C atoms is in between the 
mass of F and H atoms. 

TABLE I. List of symmetries and phonon frequencies at 
different ^points for graphane and fluorographene. R and I 
stand for Raman and infrared active modes, respectively. 



Symmetry 


Graphane 


Fluorographene 


(activity) 


r 


M 


K 


r 


M 


K 


A 2u (R) 





407 


474 





211 


260 


E„ (I) 





511 


474 





260 


260 


E„ (I) 





605 


806 





291 


298 


E 9 (R) 


1114 


1057 


1008 


250 


318 


356 


E 9 (R) 


1114 


1088 


1008 


250 


459 


510 


E u (I) 


1170 


1116 


1117 


294 


580 


510 


E u (I) 


1170 


1172 


1180 


294 


633 


961 


Ai 9 (R) 


1170 


1299 


1290 


686 


1039 


962 


E 9 (R) 


1339 


1302 


1290 


1244 


1106 


962 


E s (R) 


1339 


1317 


1322 


1244 


1116 


1114 


A 2u (I) 


2739 


2733 


2739 


1219 


1117 


1114 


Ai 9 (R) 


2783 


2746 


2739 


1312 


1188 


1150 



The 12 phonon modes of the chair conformation of 
graphane and graphene fluoride, which has the P-3ml 
symmetry (space group 164), belong to 4 irreducible 
representations. The frequency of the different phonon 
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FIG. 3. (Color online) The constant-volume specific heat of 
graphane (black solid line) and fluorographene (red dashed 
line). At high temperature, the specific heat of graphane is 
much larger than the one of fluorographene. 



modes at three high-symmetry points of the Brillouin 
Zone (r, K, and M) are summarized in Table H] E ff and 
Ai 9 are Raman active modes while E u and A2 U are in- 
frared active modes. Our theoretical values of the Ai ff 
and the two-fold degenerate E g modes of graphene flu- 
oride at the r point, 1244 cm -1 and 1312 cm -1 are in 
good agreement with the two Raman active modes found 
experimentally in Ref. H at 1270 cm" 1 and 1345 cm -1 , 
respectively. The same agreement is observed between 
the infrared active mode A2 U at 1219 cm -1 and the re- 
port of an infrared active mode at 1204 cm -1 in Ref. d. 

Using the phonon DOS we can also calculate the 
constant- volume specific heat of both materials. This is 
achieved using the expressions described in Ref. [2(| One 
should keep in mind that we work within the harmonic 
approximation, so the differences with experimental val- 
ues will become larger for higher temperatures. The hy- 
drogenation of graphene was found to be reversible at 700 
K in Ref. S while the fluorination of graphene could be re- 
verted at 870 K in Ref. 0. The calculated values of C„ are 
shown in Fig. [3] where we only report the behavior in the 
region of stability of each compound. The specific heat 
of graphane is larger than the one for fluorographene, as 
can be expected considering that the classical limit of 
the specific heat is given by the Dulong-Petit law and is 
equal to 3R/M, where R is the ideal gas constant and 
M the molar mass. Because the molar mass of fluoro- 
graphene is larger than that of graphane, the value of 
the specific heat will be smaller at high temperature. To 
our knowledge there are no experimental measurements 
yet to compare with. The specific heat of graphene at 
high-temperature [see Fig. 24 of Ref. is in between 
the specific heat of graphene fluoride and graphane as 
expected from the Dulong-Petit law. 
In summary, the experimental determination of the 
phonon dispersion relations of graphane and fluoro- 
graphene can be very useful in the characterization of 
these materials. It is worth noting, that both graphene 
derivatives are wide band gap materials. This allows 



to clearly discriminate, for example, the Raman activ- 
ity originated in fully-fluorinated and fully-hydrogenated 
regions during the synthesis of both compounds by using 
a laser with appropriate energies. 
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